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We have developed a relativistic point-coupling model of nuclear many-body dy- 
namics constrained by the low-energy sector of QCD. The effective Lagrangian is 
characterized by density-dependent coupling strengths determined by chiral one- 
and two-pion exchange (with single and double delta isobar excitations) and by 
large isoscalar background fields that arise through changes of the quark conden- 
sate and the quark density at finite baryon density. The model has been tested 
in the analysis of nuclear ground-state properties along different isotope chains of 
medium and heavy nuclei. The agreement with experimental data is comparable 
with purely phenomenological predictions. The built-in QCD constraints and the 
explicit treatment of pion exchange restrict the freedom in adjusting parameters 
and functional forms of density-dependent couplings. It is shown that chiral pionic 
fluctuations play an important role for nuclear binding and saturation mechanism, 
whereas background fields of about equal magnitude and opposite sign generate 
the effective spin-orbit potential in nuclei. 

1. Introduction 

In this work we would like to investigate the connection between QCD, its 
symmetry breaking patterns, and the nuclear many-body problem. Usual 
nuclear structure approaches 1 are consistent with the symmetries of QCD 
(in particular chiral symmetry) but only functional forms of the interaction 
terms can be determined. Model parameters cannot be constrained at the 
level of accuracy required for a quantitative analysis of structure data; they 
can only be estimated with the Naive Dimensional Analysis 2 . 
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The approach we propose is based on the following ingredients: 

(1) The presence of large isoscalar background fields which have their 
origin in the in-medium changes of the scalar quark condensate and 
of the quark density. 

(2) The nucleon-nucleon interaction is described by one- and two-pion 
exchange (with medium insertions and delta isobar excitations), in 
combination with Pauli blocking effects. 

The first point has a clear connection with QCD sum rules at finite den- 
sity 3 , in which large nucleon self-energies naturally arise in the presence of 
a filled Fermi sea of nucleons. 

The second point is motivated by the observation that, at the nuclear mat- 
ter level, the nucleon Fermi momentum kf, the pion mass m w and the A~N 
mass difference represent comparable scales 4 . Pionic degrees of freedom 
are threfore included explicitly (through density-dependent couplings), in 
contrast to the phcnomenological relativistic models, in which effects of it- 
erated one-pion and two-pion exchange are treated implicitly through an 
effective scalar field 5 . 

2. The model 

The model is defined by the Lagrangian density 3, 

£ = £free ~\~ &4f ~\~ £der ~\~ t-'coul (1) 

with the four terms specified as follows: 

Cfree = ^[ij^ - M N }^ (2) 

Us = -\g s {p)W>)W>) - \gv(p){^1^){^) (3) 
-\G T s{p){i>Ttl>) ■ {4>Tip) - ^G T v{p){4>t-i^) ■ (^tW) ( 4 ) 

Cder = -\d s {&WWW>) - \D v {d v ^^){d v Wi>) (5) 
—Dts^^I)) ■ - ^D TV (d^Tj^) • (9>T 7 ^)(6) 

C coul = -eA^Hp-i^ - \f^ v . (7) 

This Lagrangian is understood to be formally used in the mean-field approx- 
imation, with fluctuations encoded in density-dependent couplings Gi{p). 



a more details about density-dependent hadron field theory can be found in 



February 7, 2008 18:35 Proceedings Trim Size: 9in x 6in 



finclli'cortona 



Their functional dependence will be determined from finite-density QCD 
sum rules and in-medium chiral perturbation theory 

G i (p) = G^(p) + G^(p), (8) 

where the index i labels all isospin-Lorcntz structures of Eq. (1). The 
variation of the Lagrangian with respect to ip, leads to the single- nucleon 
Dirac equation 

[ 7 „(i^ - V - V£) - (M + S)]i/> = 0. (9) 

In addition to the usual self-energies V* and S, the density-dependence of 
the vertex functions produces the rearrangement contribution 13 

(IdGs o ldG TS ^ ^ 1 8G V .„ 1 8Gtv ^ 

(10) 

The inclusion of this additional term ensures energy-momentum conserva- 
tion and thermodynamical consistency. For a complete treatment of the 
density dependent point-coupling model, the reader is referred to 7 . 



3. Interactions 

3.1. Self-energies from QCD sum rules 

In leading order, which should be valid below and around saturated nuclear 
matter, the condensate part of the scalar self-energy 

is expressed in terms of nucleon sigma-mass term (ctjv = {N\m q qq\ N)) and 
the quark masses. At the same order, the time-component of the vector 
self-energy reads 




where the quark baryon density is related to that of the nucleons by 
(q' f q)p = \p- In both cases, ~ 1 GeV is the characteristic scale 
(the Borcl mass) , which approximately separates the perturbative and non- 
perturbative energy domains. 



b the four-velocity n M is defined as (1 — v 2 ) 1 / 2 (l,v) (v = in the rest-frame of the 
nuclear system) 
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For typical values of the nucleon sigma mass term ctjv (—50 MeV 8 ) 
and m u + nrid (— 12 MeV at a renormalization scale of 1 GeV 9 ), the 
in-medium QCD sum rules predict large scalar and vector self energies of 
about equal magnitude (~ 300 — 400 MeV in agreement with relativistic 
phcnomcnological models 1 ), and opposite in sign. Neglecting corrections 
from higher order condensates, these estimates have large uncertainties and 
the error in the ratio E^/Ey^ ~ — 1 is about 20%. 

Given the self-energies arising from the condensate background, the cor- 
responding equivalent point-coupling strenghts Ggy are simply determined 

by 

v (0) v (o) 

G < > = 5s_ and G^ = ^. (13) 
5 Ps V P 

At leading order the condensate terms of the couplings are constant and do 
not contribute to the rearrangement self-energy. 



3.2. Self-energies from in-medium chiral perturbation 
theory 

In recent years the nuclear matter problem has been extensively studied in 
the framework of in-medium chiral perturbation theory. The calculations 
have been performed to three-loop order in the energy density and include 
one-pion exchange Fock term, one-pion iterated exchange and irreducible 
two-pion exchange terms with medium insertions and delta isobar excita- 
tions 10 . By adjusting the coupling constants of few NN contact terms, and 
encoding short-range effects not resolved at relevant scales, to the proper- 
ties of the empirical saturation point of isospin-symmetric nuclear matter 
(Eq = — 16 MeV and po = 0.16 fm~ ), several aspects of the problem have 
been successfully investigated: the symmetric and asymmetric nuclear mat- 
ter equation of state (EOS), single-particle properties, the low-temperature 
behaviour, and the connection with non-relativistic nuclear energy density 
functionals 10,11 . 

The resulting nucleon self-energies are expressed as expansions in powers 
of the Fermi momentum kf. The expansion coefficients are functions of 
kf/m^, the dimensionless ratio of the two relevant small scales: 

^HPT = f[kf/mv) 3_ where f{kf/mv) = J2c J*L) . (14) 

1V1 N k V m 7T/ 

The density-dependence of the strength parameters is determined by equat- 
ing the point-coupling self-energies in the single-nucleon Dirac equation (9) 
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with those calculated using in- medium chiral perturbation theory : 

G P PS = Zr PT (k f ) (15) 
G^p + V^=^ HPT (k f ) (16) 

G^ipi = ^r T (k f ) (it) 
4tp 3 = ^r T (k f ). (is) 

In Fig. (1) the resulting equation of state of isospin-symmetric nuclear mat- 
ter is compared with the CHPT nuclear matter EOS 10 . The agreement is 
satisfactory, and the small difference can be attributed to the approxima- 
tions involved (i.e. the momentum dependence is neglected). 




Figure 1. Binding energies for symmetric nuclear matter as a function of baryon density. 
The solid curve (CHPT) is the EOS calculated by using in-medium CHPT. The EOS 
displayed by the dashed curve (CHPT-mapping) is obtained when the CHPT nuclcon 
self-energies are mapped on the self-energies of the relativistic point-coupling model with 
density-dependent couplings. The dotted curve denotes the latter EOS when rearrange- 
ment terms arc neglected. 



4. Results 

In this section the relativistic point-coupling model will be applied to the de- 
scription of ground-state properties of finite nuclei. The density-dependent 
coupling strenghts will be determined in a least-squares fit to observables 



Tor the definition of p, p 3 , P3, p s 3 the reader is referred to 



7, 2008 18:35 Proceedings Trim Size: 9in x 6in 



finclli'cortona 



6 



(O) of a rather large set of nuclei along the valley of /^-stability (from lb O 
to 214 Pb): 



X 



E E 

nuclei \observ. 



{Oth — Oexp) 



weights 



(19) 



In Fig. 2 the overall agreement is shown in comparison with two stan- 
dard relativistic meson-exchange parametrization: NL3 12 , with explicit 
higher order self- interaction terms for the cr-meson, and DD-MEl 13 , with 
phenomenological density-dependent couplings strenghts. The interaction 
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Figure 2. Relative errors for the observables binding energies (upper panel), and charge 
radii (lower panel), in comparison with the results of NL3 and DD-MEl. 



determined by the least-squares fit procedure contains pionic fluctuations 
on top of the large background fields. Nonetheless, it is interesting to an- 
alyze the role of pionic fluctuations, without the presence of background 
fields, and how these fields modify the single-particle spectra in finite sys- 
tems. In Fig. 3 we display the single proton levels in 40 Ca without and with 
background fields. It is obvious that the the spin-orbit potential 

V °^^{l^ V °^- S ^) l - S > (20) 

represents a short-range effect (see also 10 ), determined by the condensate 
structure of QCD. The isospin-dependent part of the interaction has been 
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Figure 3. Single proton levels for 40 Ca with only pionic fluctuations (left panel) and 
with the inclusion of background fields (right panel). In the latter case the degeneracy 
of the p- and d-levels is removed. 



studied in the description of the neutron radii. Standard rclativistic mean- 
field calculations sistematically overestimate the difference between neutron 
and proton radii 1 . It turns out that the density dependence of the isovector 
channel of the interaction is crucial in order to reproduce these observables, 
as also shown in 13 . In Fig. 4 we show the calculated values of r n — r p for the 




Figure 4. Predictions (diamonds) for the differences between neutron and proton radii 
of Sn isotopes, in comparison with NL3 (empty circles), DD-ME1 (empty triangles) 
predictions and experimental data (empty squares) . 



Sn isotope chain in comparison with experimental data 1 . The agreement 
is excellent, and it is a clear indication that the isovector channel can be 
successfully described by pionic fluctuations. 
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5. Conclusions 

It has been demonstrated that an approach to nuclear dynamics constrained 
by the low-energy sector of QCD provides a quantitative description of 
properties of nuclear matter and finite nuclei. 
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